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LHC hints could come from the Higgs sector, where loop corrections involving (j) lead to 
significantly modified /i — > 77 and h — > 7Z branching ratios. 
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1. INTRODUCTION AND MOTIVATION 

Recent analyses of Fermi-LAT data have revealed a line-like feature in the cosmic gamma ray 
energy spectrum from the Galactic Center at an energy ~ 130 GeV [IHl]. Additional hints for 
a 130 GeV photon line were seen in galaxy clusters {5} and unassociated Fermi-LAT sources [B] 
(see, however, [ZHS])- At present, it is not clear whether these features are due to an instrumental 
effect or due to physics beyond the Standard Model (SM). Validation tests done in the original 
Refs. [IHS]; as well as additional checks using the public data performed in Refs. [TOHn] . have so 
far not identified an obvious problem with the data, but an official analysis by the Fermi-LAT 
collaboration will certainly shed further light on the issue. 

In this paper we assume that the signal is evidence for dark matter (DM) particles x annihilating 
into two photons, XX ~^ 77) oi a photon and a Z boson, XX ~^ -^7- In the former case, the DM 
would need to have a mass ~ 130 GeV and an annihilation cross section {cr{xx — ^ 77)'yrei) ~ 
1.3 X 10~^^cm^/s [2], whereas if the signal is due to the annihilation process xX ~^ -^7) one obtains 
My, ~ 144 GeV and (crurei) ~ 3.1 x 10~^''cm^/s jlOj . The fact that we see a photon signal requires 
that DM couples to a state (p that is charged under the electroweak gauge group. Annihilation can 
then proceed through (j) loops. The required annihilation cross section (XWrcl is roughly an order of 
magnitude smaller than what is required for a thermal relic, but still large for a loop suppressed 
process. It is, for instance, much bigger than what is expected from a singly charged particle (j)'^ 
running in the loop, unless the coupling of (p'^ to DM is large, close to the perturbativity limit |15] . 
Additionally, if the charged particles <j)^ in the loop are lighter than rriy., the DM can annihilate 
into them at tree level. These annihilations would contribute significantly to the continuum photon 
emission from the galactic center due to final state radiation and decays of secondary pions. The 
resulting annihilation cross sections are typically excluded by strong bounds on the continuum 
photon emission from the galactic center |16H20j . 

Many models have been proposed to circumvent these problems |20H58j . In this paper we focus 
on a particular set of models that can lead to interesting signals at the LHC. In these "secluded 
dark matter" models, DM couples to the visible sector primarily through loops of a new electroweak 
multiplet (j). For concreteness we focus on examples where i;^ is a scalar with vanishing vacuum 
expectation value (vev). The salient features of this type of model are 

• The DM annihilation cross section to photons is enhanced because some states in the medi- 
ator multiplet carry large electric charges. 

• For suppressed DM-Higgs coupling, the continuum photon bounds are avoided because then 
the dominant annihilation to W and Z bosons is generated at one loop, and to SM fermions 
only at two loops. The correct relic density is obtained if (j) is somewhat heavier than the 
DM. 

• If the mediator (p couples to the Higgs boson h, the branching ratios for the decays /i — )• 77 
and h — )• 7Z are altered. If the new particle discovered recently by the ATLAS and CMS 
collaborations [HUl [50] is indeed a SM-like Higgs boson, the experiments could see these 
modified branching ratios in future precision measurements. 

• The charged components cp"^^ of the mediator multiplet cp are produced at the LHC through 
their electroweak gauge couplings. Decays to off-shell W bosons lead to multi-lepton final 
states, however, for large parts of the parameter space the leptons are so soft that the signal is 
not observable at the LHC or the Tevatron. The best probe of (p"^^ production are then final 
states with a photon and large missing energy (plus possibly other visible particles) because 
of the large couplings of (p^^ to the photon. Moreover, for very small mass splittings, the 
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lifetimes of the cp are so long that they can appear as anomalous charged tracks in the 
inner detectors of ATLAS and CMS. 

• All 4-scalar couplings are perturbative and continue to be so up to the Planck scale. In 
particular, the DM-(j) coupling can be relatively small and still lead to a large gamma ray 
signal because of the large c/)"^ charges. For the large SU{2) representations {N > 3) 
considered here, the weak gauge coupling becomes non-perturbative below the Planck scale, 
see e.g. |61j . This implies that perturbative grand unification is only possible if the model is 
embedded into a more complete theory at an intermediate scale. For N < 9, the embedding 
(or, alternatively, non-peturbativity of the weak interaction) does not have to occur at scales 
below several 100 TeV, outside the reach of the LHC. 

The connection between the 130 GeV gamma ray line and an enhanced /i — )• 77 signal at the 
LHC has been made also in [36] for a model with an electroweak triplet mediator. While in |36j 
implications for other LHC searches were not elaborated on, we keep the discussion as general as 
possible and explore also LHC signals aside from the enhanced Higgs to diphoton rate. We also 
consider general electroweak multiplets beyond the triplet, but for numerical examples we will use 
electroweak quintuplets as mediators. We will discuss to what extent electroweak multiplets are 
constrained by precision Higgs physics, by searches for anomalous charged tracks, and by monojet, 
monophoton and photon + MET + X searches. In the context of the 130 GeV gamma ray line, 
LHC final states with a photon and missing energy were also considered in [55] in the context of 
models with Z' and axion mediators. Since in these models, the photon is produced as part of 
the hard process, mono-photon searches are more constraining than in our models, where photons 
are only produced radiatively. Finally, independently of the 130 GeV line, the effects of a scalar 
electroweak quartet on Higgs boson decays to 77 and Zj have been considered previously in [62j. 

The paper is organized as follows. In section [2] we introduce the class of models we consider in 
more detail. Section |3] focuses on DM annihilation into photons and on continuum photon emission 
bounds. In section |4] we discuss the cosmological history and prospects for DM direct detection. 
Section [5] deals with existing electroweak precision constraints. The collider phenomenology is 
discussed in section [6], and the modifications of the Higgs boson properties in section [7] We 
conclude in section |8j while calculational details are relegated to the appendices. 

2. MODEL SETUP 

We consider an extension of the SM by a scalar A^-dimensional SU{2)l multiplet (j) of hyper- 
charge Y^. If > 5 there are no renormalizable couplings to the SM linear in (f)} The (j) fields 
then interact with the SM only through Higgs portal and gauge interactions, 

C D \D^cl)\^ - m\^^ - X^h<PUh^H - X'^Hi<P^T^<P){H^T^H) - X^{<pUf, (1) 

where and are the generators of the SU{2)l representations N and 2, respectively (their 
normalization is given in Appendix [Aj) . We assume the Higgs portal coupling X^f^u to be either 
positive or negative but with |A(^hz;^<C m^, so that cj) does not develop a vacuum expectation 
value. Here, v is the vacuum expectation value (vev) of the Higgs. For the same reason (and other 
reasons discussed below), X'.^ should not be too large in magnitude. Expanding the covariant 



^ This is also true in the case that A'^ < 4, except for some specific values of Y^. 
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derivative in eq. ([T]) gives interactions between (j) and the electroweak gauge fields, 

Note that since eq. ([T]) is the most general renormahzable Lagrangian, the Z2 symmetry (f> ^ —cj) 
is accidental. The neutral component of (j) can thus be stable and a DM candidate in principle. 
However, the annihilation process — )• W~^W~ , which occurs at tree level, has too large a 
cross-section to give the observed DM abundance today for (p masses below TeV scales [631 El] • 
Moreover, even with the correct relic density, the same annihilation process for relatively light 
(f) in the present day would be in tension with observations of dwarf galaxies by the Fermi-LAT 
collaboration [65]. In order to accommodate the tentative Fermi-LAT line, we therefore introduce 
an additional real vev-less SM-singlet scalar x, which has direct couplings only to the other scalars: 

C D Id^xdf'x - \mlx^ - \hX^H^H - A^^xW- (3) 

Here, we have introduced by hand a Z2 symmetry that stabilizes Xi assumed that is 
chosen such that x does not develop a vev. If we wish to explain the tentative Fermi-LAT gamma 
ray line at ~ 130 GeV, the DM mass is fixed at ~ 130 GeV or ~ 144 GeV, depending on 
whether decays to 77 or 7Z are dominant as we discuss below. In general, however, can be 
arbitrary. The mass parameter of the weak multiplet, m^, is also free, but as we will see below, 
phenomenologically most interesting is the region where > M^. After electroweak symmetry 
breaking, the physical masses of the particles in the (j) multiplet, M^n±, receive an additional 
contribution from the Higgs vev, v = 246 GeV. The coupling in ([T| leads to an overall shift, 
while the X'^^ term after EWSB, 

- A;^(</.tr- 0)(/7t^a^) ^ + 1 A'^^^2^tr3 <^ , (4) 

leads to a mass splitting 

A(m2) = --)^^Hhv\ (5) 



between the (j) component with eigenvalue and the = component. There are three inter- 
esting regimes of the \'^^ coupling; for \'^^ ~ the splitting is tens of GeV, for \'^^ ~ 0(0.1) 
the splitting is several GeV, while for \'^^ = a splitting arises only from one loop electroweak 
corrections and is tens to hundreds of MeV as shown in fig. [T] (right). For large mass splittings 
the decays of charged particles are easily observable at the LHC and are excluded, so we will 
be interested in smaller values of of 0(0.1) or below. The parameter in general cannot 
be made arbitrarily small without fine-tuning since it can be generated from a loop with A"^ and 
B^ on the two internal lines. For = 0, however, this contribution is zero (cf. eq. Q) so that 
\'^^ ~ is natural in this case. For 1^ 7^ there is a log divergent contribution to the bare X'^^j 
coupling. Even if the cut-off of the theory is at the Planck mass, however, such a contribution is 
only log(Mp;/Mvi/)2Y'0a/47r ~ 0.1, so that the values of A'^^ chosen in fig. [T] (left) are natural. 

As already mentioned, for small values of A^^ an important contribution to the mass splitting 
are the 1-loop electroweak radiative corrections. The resulting mass splitting is given by [63j 

Mq - Mq, = ^^{(Q^ - Q")sl,f{Mz/M,) 

+ {Q- Q'){Q + Q' + 2Y^)[f{Mw/M^) - f{Mz/M^)]], 
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Figure 1: The mass difference in MeV between the singly charged and neutral mediators, 
function of the 0" mass M^o and the coupling constant X'^j^ defined in eq. (fll) 
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where Q and Q' are the electromagnetic charges of two component fields of (p, a is the fine structure 
constant, is the hypercharge of the multiplet, and sw = sin 6\y is the sine of the Weinberg angle. 
The loop function / is given by 



fir) 



2r^ In r + (r^ - 4)^/^ In a] /4, with A = (r^ - 2 



4)/2. 



(7) 



Here, the UV divergence has been absorbed into the renormalization of and X'^jj- (we are using 
the scheme A; = in the notation of [63j). Numerically, /(I) = 2.72, so that for ~ 0(100 GeV) 
the mass splitting due to electroweak corrections is tens of MeV. 

Because of the accidental Z2 symmetry in eq. ([T]), the lightest component of (p is stable. If 
we view the model only as a low energy effective theory, however, the Lagrangian in eq. ([T]) is 
supplemented by higher dimensional operators which can allow the lightest component of (j) to 
decay. If (p forms an A^-dimensional multiplet, then the lowest dimensional operator mediating this 
decay needs to contain at least — 1 SM doublets. For example, the choice = 5, V^ = 2 allows 
us to include the operator 



(8) 



For general N, operators of this type will be suppressed by 1/A^~^, where A is the cut-off scale of 
the effective theory. 

In the following, we consider in detail two benchmark cases, one in which we assume that the 
lightest component of the (f) multiplet is stable on cosmological timescales, and one in which it 
decays rapidly through higher dimensional operators, see table |lj In the stable case, the lightest 
component of (j) contributes to the dark matter relic density at the subdominant level. For instance, 
for the benchmark model listed in the left part of table |l| its relic density is Vt^oh? = 3.6 x 10~^. 
It is thus important that the lightest component of (j) is electrically neutral and does not couple 
to the Z — if it did, its scattering cross section on nuclei would be in conflict with direct detection 
constraints. To avoid couplings to the Z, we have to ensure that (fP has = 0, which is only 
possible for odd multiplet order A'" and requires Ys = 0. We choose A^ = 5 for definiteness, but 
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Table I: The input parameters, resulting mass spectra and relic densities for the two benchmark points: 
A — 5-plet with stable cj/' (left), and a = 2 5-plet with cj/' allowed to decay through higher- 
dimensional operators (right). Note that each 0"* state is associated with an antiparticle <^*'"* carrying 
equal but opposite charge. 

larger multiplets are also viable. To make sure that (p^ is indeed the lightest component of cp, we 
also assume that X'^^j- is small enough that the mass splittings among the components of (p are 
dominated by electroweak corrections. For YJ^ = these lead to small positive mass shifts for the 

7^ charged components compared to the neutral one. We also set Xy^n = so that there is 
no XX h ^ WW annihilation at tree level. The phenomenological consequences of relaxing this 
assumption will be addressed below 

If (p can decay through higher-dimensional operators, there are much fewer constraints. For 
example, if the decay is fast enough, all components of (p could be charged and there is no con- 
straint on which component is the lightest one. Here, we will nevertheless assume that the lightest 
component is electrically neutral. The complete set of model parameters for the two benchmark 
cases is given in table |lj In both of them, we focus on = 5 multiplets, but we will also comment 
on higher multiplets below. 

3. GAMMA-RAY ANNIHILATION SIGNAL 

We are now ready to discuss in detail the phenomenology of the DM models introduced in 
section [2| where DM-SM interactions are mediated by a scalar SU(2) N-plet. We begin by con- 
sidering indirect detection constraints, in particular the signals of DM annihilation in the gamma 
ray sky. As mentioned in the introduction, one of our motivations is the tentative line-like feature 
observed in Fermi-LAT gamma ray data from the galactic center and other DM-rich regions in the 
sky [IHS]. This signal, as well as possible gamma ray lines that may be discovered in the future, 
can be due to either xx ~^ 77 or xx ~^ 1^ annihilation. The process XX ~^ 7^ is not generated 
in the models we consider because the initial and final states would have different C parity. Both 
XX — ^ 77 or XX ~^ 1^ proceed through diagrams of the form shown in fig. [2] In both of our 
benchmark points from table [l] we have X^h = so that only the topologies fig. M (a) and (b) 
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Figure 2: Representative diagrams contributing to xx ^ Hi 7^1 ZZ annihilations 



contribute. The annihilation cross sections are then given by 
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with /3 = M^/M^, 7 = 4M^/M|, and the loop functions 
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for X > 1, 
for X < 1. 



for X > 1, 



yjx — 1 arcsin \/x ^ for x < 1 . 



(11) 



(12) 



In the above expressions, sw and cw denote the sine and the cosine of the Weinberg angle, respec- 
tively, Q and Is are the electric charge and the third component of the weak isospin of the (j) held 
components, and the sums run over these components. 

The predicted values of the annihilation cross sections (crfrei)^^ and (ciUrel)^^ are shown in 
fig. 3 These should be compared to (cjfrci)^^ = (l-27!to'43) x lO"^^ cw? /s and (crurei) 



(3.14I^;99) X 10^27 cm^/s^ which, for the Einasto DM profile were shown in [21 EO] to explain the 
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Figure 3: Contours of constant annihilation cross section (crfici) for the annihilation processes XX ~^ 11 
(blue solid lines) and xX ~^ 1^ (red dashed lines) as a function of M^a (the mass of the neutral component 
of the mediator multiplet (f) and A^^^ (the coupling of DM to 0) . Motivated by the two benchmark models 
given in table |lj we take 4> to be an SU (2) 5-plet with = and no isospin-dependent couplings to the 
Higgs (A'^^ = 0) in the left panel, whereas in the right panel we chose = 2 and X'^^j = —0.1. Our choice 
of DM mass, = 144 GeV for xX ~^ 1^ and = 130 GeV for xX ~^ 11 is motivated by the tentative 
Fermi-LAT gamma ray line signal [IH31 [S] . The thick red line denotes the values of M^o and A;^^ for which 
the correct DM relic density SIdm^^ = 0.112 W] is obtained if all the other model parameters are fixed as 
m table in The error on the relic density from WMAP, ±0.0056, is below the resolution of the plot. 



Fermi-LAT feature at 130 GeV for DM masses = 130 GeV and = 144 GeV, respectively. 
We see that annihilation cross sections > 10^^^ cm^/s are easily obtained in our model for A^,^ well 
within the perturbative regime, and without the need for tuning between M^o and M^. Notice 
the qualitative change in the dependence of {(yv^f,\)^^ on M^o and A^,^ when M^o approaches M^. 
The reason is that for M^- < M^, the loop diagrams in fig. [2] acquire an imaginary part because 
direct annihilation xX ~^ becomes possible. This effect is much more pronounced for the stable 
benchmark point (fig. [3| left plot) due to the near-degeneracy of the components of (j). In the 
unstable case, the non-zero hypercharge assignment allows for destructive interference in xX ~^ ^1 
for M^o ~ 125 GeV, leading to the "kink" visible in the right plot of fig. [3j 

To illustrate the dependence of (cTVrei).y^ and {(TVrei)jz °^ quantum numbers of (p, we show in 
fig.jijleft (right) contours of constant (ct'rcl)-^^ (-y^) = 3.14 (1.27) x 10^^'^ cm'^/s, i.e. for the central 
values of annihilation cross sections motivated by the tentative gamma ray line at 130 GeV [2l 110]. 
Several different choices for the multiplet dimension N and its hypercharge are shown. As 
expected, it is easiest to obtain the annihilation cross sections required to explain the 130 GeV 
line in models with large N and thus highly charged component fields of (p. Note that {crvrei}^^ 
increases with because higher charge states appear for large F^, whereas (crfrei)^^ decreases 
with because of stronger cancellation between I3 and s^^Q in the term at the end of the first 
line of eq. (10). 

Besides the annihilation to two photons, the DM in our model also annihilates to W~^W~ and 
ZZ. If we set the DM-Higgs coupling X^h to zero, as in our benchmark points from table |l| 
annihilations to W~^W~ and ZZ first occurs at 1 loop level. The annihilation cross section is 
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Figure 4: Contours of constant annihilation cross section {<yvic\)^z ~ i-^^ x 10^^^ cm'^/s (left) and 
(cruici)^^ — 1.27 X 10^^^ cm'^/s (right), motivated by the tentative Fermi-LAT gamma ray hue [IHSl [5], 
as a ftmction of M^o (the mass of the neutral component of the mediator SU{2) multiplet (jj) and A^^ (the 
coupling of DM to 0). The results are shown for multiplet sizes = 3, 5, 7, 9 (green, orange, magenta and 
blue lines), and hypercharge = 0, 1, 2 (solid, dashed, dotted lines). The choices of DM mass, = 144 for 
XX ~^ iZ (left panel) and — 130 GeV for XX ~^ 11 (right panel) are also motivated by the Fermi-LAT 
line. The remaining input parameters are as in the corresponding columns of table |lj The yellow bands for 
the benchmark models A = 5, 1^ = (left) and A = 5, y0 = 2 (right) show the la experimental ranges 
for (CTWroi)^z = (3.14t°||^) X 10-27 cm^/s and (crwrci)^^ = [1.27^0%) x lO'^'^ cmVs, respectively (obtained 
using an Einasto halo profile in [21 [IQ]). The thick red line denotes the values of M^o and A^^ for which 
the correct DM relic density Oqm^^ = 0.112 69 is obtained. The error on the relic density from WMAP, 
±0.0056, is below the resolution of the plot. 

then smaller than the bounds from continuum gamma rays in Fermi-LAT. Using FeynArts |70], we 
estimate that for the benchmark point with the stable 5-plet (left part of table |l]), the annihilation 
cross section to W^W~ is (aurci) = 2.0 x 10~^^ cm'^/s, and the one to ZZ is (crvrci) = 5.2 x 
IQ-^^ cm^/s. In the case of the unstable 5-plet benchmark point (right part of table |l]), the 
annihilation cross section to is {crv^ex) = 5.3 x 10^^^ ctd? /s and the one to ZZ is (cri'rel) = 

2.6 x 10^^'' cw? /s. The bound from continuum gamma rays from the galactic center is (ufrei) = 

2.7 X 10-2^ cm^/s for annihilation to W'^W and (avrci) = 3.2 x 10-2*^ cm^/s for the ZZ final 
state J9] . The continuum photon constraints can also be translated into a constraint on which 
is X^H < 0.03. 

4. RELIC DENSITY AND DIRECT DETECTION 

We now investigate the dynamics of DM freeze-out in the early Universe for the class of models 
given by the Lagrangians ^ and At very high temperatures, the DM x is kept in thermal 
equilibrium through two channels: (i) s-channel Higgs exchange xX ^ h ^ WW, ZZ and (ii) 
direct coupling to the mediator field (/>, xX ^ 'P4'- 4': i^i turn, is kept in thermal equilibrium with the 
SM particles through its electroweak interactions. The amplitude for process (i) is proportional 
to the coupling constant XxH^ which is constrained by the requirement that secondary gamma 
rays from DM annihilations in the Galactic Center today should not overshoot the Fermi-LAT 
constraints on the gamma ray continuum. Since generating the correct DM relic density Qh'^ = 
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0.1120 lb 0.0056 |69j through XX -H- /i o WW,ZZ alone is only marginally allowed, we will not 
entertain this possibility here. Instead, we focus on the case where the correct relic density of DM is 
determined by the "forbidden" annihilation channels [3 H 1721 [73] , xX ~^ These channels are not 
kinematically accessible for nonrelativistic DM since < M^. Therefore, they do not contribute 
to DM annihilations today, avoiding indirect detection constraints. In the early universe, however, 
they can still be effective if is not too much larger than M^^ so that xX ~^ 4'4' is still accessible 
from the high-energy tails of the thermal DM energy distribution. Depending on the quantum 
numbers of (j), we find that the mass gap — required to explain the observed relic density 
is several tens of GeV to 100 GeV. Using MicrOMEGAs [71] we estimate that for our N = 5, 1^ = 
benchmark point (left part of table|l]), DM freeze-out occurs at M^/T ~ 25. The components of the 
multiplet (j) remain in thermal equilibrium until M^/T ~ 26 (which is equivalent to M^i^/T ~ 33) 
and have lower relic density than x- Around freeze-out the DM velocity is ~ 0.3 c so that it can 
still annihilate into the slightly heavier (p. This maintains equilibrium with the thermal bath and 
as a result, the DM freezes out at around the same time as it would in the simple thermal WIMP 
scenario, giving the correct thermal relic density. 

Prom the thick red lines in figs. [3] and |4] we see that, for both the stable and unstable 5-plet 
benchmark models, the correct relic density can be obtained for a number of different parameter 
choices. The calculations were performed using MicrOMEGAs [71], but cross-checked by solving 
the relevant Boltzmann equations numerically in Mathematica. In the case of the stable the 
neutral component i;^'^ constitutes part of the DM relic density, but as mentioned in section [2| its 
abundance is expected to be small because of its efficient annihilation to W^W~ and the resulting 
late freeze-out. Indeed, using MicrOMEGAs, we find Q.^Qh? = 3.6 x 10~^, which is three orders of 
magnitude lower that the total DM relic density. 

We next discuss direct detection signals in the models we are focusing on. DM interactions with 
nuclei are described by the following effective operators: 

= %X'F,.F^'^ + ^ ^,M,x'm. (13) 

4> g <f> 

where in the interactions with quarks we have already included the required quark mass suppression 
due to a chirality flip. Por the suppression scale, we have chosen the (j) mass for later convenience, 
while Cp and Cq are dimensionless Wilson coefficients. The effective cross section per nucleon for 
DM scattering on a nucleus {A, Z) is 



- ' ' {-Z'Qom^,+A^m,f^-^y, (14) 



A, 



where A, Z are the atomic mass and atomic number of the nucleus, respectively, the nuclear coher- 
ence scale is Qo = \/6(0.3 + 0.89^1/3)-! fm"^ [29], and we neglected the momentum dependence 
in the electromagnetic form factor J^(|q|), replacing it with -F(0) = 2/^/tt. The average for the 
matrix elements /^'^ of the scalar operators in the second term is over neutrons and protons in the 
nucleus, where we use the values given in |75j and take proton and neutron masses equal. In our 
model Cp arises at 1-loop from (j) running in the loop, while Cq arise at two loops involving (p and 
Z, W exchanges. Numerically, the typical size of the scattering cross section on Xe is 

af/ -59. 10-^« cm^ f 160GeV\ ' / 140GeV\ ' ( Cq \ ' 
a^_5.9 10 cm J V 1/16^2 +^-^1/(16^2)2; ' ^^^) 

where for simplicity we have assumed that Cq is independent of quark flavor. This is not entirely 
correct in our models, where we have Ct ~ C„ c = 1-9 x 10~^ and ^ & = 2.1 x 10~^ for stable and 
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unstable 5-plet benchmark points, respectively. (Here we have evaluated the two loop integrals in 
the limit ^ M\y^z-) For the diphoton operator the numerical values of the Wilson coefficients 
are Cp ~ 10~^. The analytical expressions are given in Appendix[Bj The numerical values should be 
compared with the present XENONIOO bound, which for ~ 140 GeV DM is afj < 4-10-^^ cm^ [75] . 

Note that the Cq Wilson coefficient also receives a nonzero tree level contribution due to single 
Higgs exchange, C^^'^^ = —X^nM^/Mfj. The coupling X^h is bounded from the continuum photon 
flux to be X^H ^ 0.03 (see Appendix [D]) , and from direct detection it is bounded to be also below 
X^H ^ 0.03. In our benchmark points we set X^h = 0. 

5. PRECISION ELECTROWEAK CONSTRAINTS 

With the addition of a new charged multiplet at scales not far above the electroweak-breaking 
scale, we might worry that there will be severe constraints from current experimental bounds on 
precision electroweak observables. Because our new multiplet does not couple directly to any SM 
fermions or induce any new symmetry breaking, we expect no significant corrections to flavor 
physics observables or processes such as anomalous electric dipole moments. We will therefore 
restrict our attention to observables linked closely to the gauge sector, namely the running of 
the gauge couplings, the S, T, and U parameters, and contributions to the anomalous magnetic 
moment {g — 2) of the electron and muon. 

The presence of new charged matter, particularly in a high representation of the gauge group, 
can have a significant impact on the running of the gauge couplings g' and g (corresponding to 
U{1)y and SU{2)i, respectively.) If the rate of running is increased substantially, the gauge 
couplings can become non-perturbative at relatively low energy scales. We will not insist on 
perturbativity of the couplings all the way up to the Planck mass, but only up to multi-TeV cutoff 
scales, ensuring that our theory is valid at LHC-accessible energy scales. 

The contribution of the multiplet (j) to the one-loop /3-function coefficient hi is given by |61| 

h^^i = ^T,(N), (16) 

where for U{\)y the quantity Tj(N) = Y^, while for SU{2)i it is equal to the trace invariant 
C(N) = (A'^^ — N)/12. For our stable benchmark case = 0, while in the unstable benchmark 
it is set to = (A^ — l)/2. In either case, the more stringent constraint on perturbativity comes 
from the SU{2)l running, due to the stronger scaling of b^^i with A^. For the choice A^ = 5, we 
have b(j)^i = 10/3, and the gauge coupling remains perturbative up to the Planck scale [61]; higher 
representations would give stronger running, with A^ = 9 leading to a breakdown in perturbativity 
at a scale on the order of 1000 TeV. 

We turn now to the "oblique parameters" S, T and U [77], which encapsulate generic contri- 
butions of new electroweak-charged physics objects to low-energy observables. These parameters 
are related to the transverse vacuum polarization of electroweak gauge bosons li ab{p)- From [77], 
eq. (3.12), the parameters are given by 

a5^4e2[n'33(0)-n'3Q(0)], (17) 

«^=^2j^[nii(o)-n33(o)], (18) 

aC/^4e2[n'n(0)-43(0)]. (19) 

Considering S first, for a scalar particle and using the definition Q = + y, the polarization 
amplitudes can be split diagram by diagram: 

s (X n'33(o) - (n'33(o) + 4^(0)) ^ -n'3y(o). (20) 
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If there is no significant mass sphtting within the multiplet, then this amphtude is proportional 
(at leading order) to tr(r^), which vanishes identically in any representation. However, in the 
presence of such a mass splitting, there will be in general a non-zero contribution to S. Computing 
the relevant one-loop amplitude, we find that 



3 sin 6,,, cos 9. 



E 

/3>0 



-^alog 



Ml 



y.t>-h. 



(21) 



where the sum is over all positive Is values for the multiplet (/>, and M^^q denotes the mass of the ^ 
state with electric charge Q. Details of the calculation are shown in appendix [C} For the unstable 
(f) 5— plet benchmark point given in table m we find the contribution S ~ —0.016. This is not large 
enough to cause any tension with the experimental constraint S = 0. 00^*^0' Jq [78], but in principle 
a larger multiplet with a substantial mass splitting could be constrained by S. 

The calculation of the T parameter is similar, but slightly more involved; details are given in 
appendix [C| We find 



—a 



^s(iV - s)5s+i,s' + \{s - 1){N - s + l)(^s-i„ 



(M2+M2) 



2Mf 



Ml - Ml 



log(M,VM2) 



(22) 



Evaluating this expression numerically for our benchmark points yields T ~ —0.013 for the stable 
case, and T w 0.0062 for the unstable case, both well within the current experimental bounds 



0.02 



+0.11 
-0.12 



1781. 



Because there is no direct coupling to standard-model fermions, contributions of the new multi- 
plet (p to the anomalous magnetic moment [gg — 2) of the charged leptons (£ = e, /i, r) will appear 
starting at two-loop order, through modifications of electroweak vacuum polarization and vertex 
functions. Because the new sector does not induce any additional breaking of gauge symmetry, in 
the limit of M^ — )• oo we expect that all such corrections must vanish due to gauge invariance. The 
leading contribution to (gi — 2) should thus scale as 1/M^. Naive dimensional analysis then gives 
us the rough estimate 



Aa^ = {Ag^ - 2)/2 



g 



Ml 



a 



Ml 



2(167r2)2 Ml 327r2 M' 



2 ' 



(23) 



which for > 100 GeV gives a contribution of about < 2 x 10^^^, three orders of magnitude 
below the current experimental uncertainty in \79\ |80] . The expected deviation of the electron 
{g — 2) from its SM value is even further from being experimentally constrained. However, the above 
estimates do not include a prefactor due to the charges and multiplicity of the (p components, which 
could easily be O(IO^) or larger depending on the exact choice of representation and hypercharge. 
Although we will not attempt it here, a precise two-loop calculation of the contribution from this 
new multiplet to would be interesting, and could potentially yield a contribution large enough 
to explain the current discrepancy between theory and experiment in this quantity |80| . 



6. COLLIDER PHENOMENOLOGY 



The only direct coupling between the dark matter x ™d the SM particles in the class of models 
discussed here is through the Higgs portal operator xxH^H, see eq. However, we have seen 
at the end of section |4] that the corresponding coupling constant should be small in order to 
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Figure 5: The LHC pair production cross section for the new electroweak multiplet (j) at our benchmark 
points. The blue and red curves are for the — Q and 1^ = 2 benchmark points, respectively. Solid curves 
are for y's — 14 TeV, dashed ones are for the 8 TeV LHC. The width of the colored bands indicates the 
theoretical uncertainty of our predictions, estimated by varying the factorization and renormalization scales 
in MadGraph by a factor of 2. 

avoid constraints from direct detection and from continuum photon emission in DM annihilation. 
Therefore, we do not expect the DM production cross section at the LHC to be large enough to 
be discovered anytime in the near future. On the other hand, the components of the mediator 
multiplet (j) can be produced abundantly at the LHC through their large electroweak couplings. 
Their decay phenomenology will depend crucially on the mass splittings between them, and since 
these mass splittings can be quite small, very interesting collider signatures are expected. We will 
now discuss the collider phenomenology of the new electroweak multiplet (f) in more detail. 



6.1. 



production and decay at the LHC 



At the LHC, the electroweak multiplet (j) would be produced mostly in Drell-Yan pair production 
processes. In fig. [5| we show the expected production cross sections at center of mass energies of 
8 TeV and 14 TeV. We see that, especially for relatively light (p {M^ < 200 GeV), as at our 
benchmark points, the production cross section is fairly large, on the order of 1 pb at y/s = 8 TeV 
and up to more than 10 pb at ^/s = 14 TeV. Nevertheless, detecting (j) is challenging because 
by assumption its lightest component is typically electrical neutral, and the decays of the heavier 
components are very soft. 

In particular, the charged components of decay via i/)"^ — )• + W* , where the off-shell 

W* gives leptons or hadrons in the final state. The relevant decay rates are |63j 



r(0"± ^ </.("-i)± + 7r±) 



(AM)2 



(24) 
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r(,...^,..-..^,e^y.). (^^-;)g(^^^)' . (25) 

In these expressions N is the dimension of the SU{2) representation of n = . . . {N — 1) /2 labels 
the components of (j), V^d is the CKM matrix element, ~ 130 MeV the pion decay constant, 
the pion mass, and AM the mass difference between and (p^"'~^^^. We have made the 



approximation that ^> AM, M^^. In eq. (25), we also set the electron mass to zero. In the 

case of the (/)"^ decay to a muon and a neutrino, (/)"^ — )• (^^""^^^/U^V^, a similar approximation, 
'Ti'/i — 0, is not appropriate, and the analytic expression for the corresponding decay rate is lengthy. 
We instead used CalcHEP [81] to compute the decay rate T numerically. 

It is clear that for small mass splittings AM the hadrons or leptons produced in (p^^ decays are 
very soft and are thus undetectable in the LHC's high energy, high luminosity environment. For 
instance, for ~ 150 GeV, even relatively large mass splittings of 5 GeV lead to a lepton with 
PT > 10 GeV in only about 2.6% (3.2%) of </) pair production events at the 8 TeV (14 TeV) LHC. 
A jet with pt > 25 GeV is produced in only 4.4% (5.9%) of the events. In these percentages, jets 
from initial or final state radiation are not included. 

If other energetic final-state particles are present, the cascade decay products can be boosted 
and therefore easier to detect. For example, requiring a final-state photon with px > 80 GeV 
leads to leptons with px > 10 GeV in 7.5% of (j) pairs produced at the 14 TeV LHC. However, 
the production cross section is reduced to 20 fb. Existing searches for e.g. W + ^ + MET |82j are 
therefore not constraining, and even future searches in this channel would be challenging, although 
a search strategy with more sophisticated kinematic cuts may be more sensitive. 

If AM is smaller than M^^ the hadronic decay modes are kinematically forbidden and only 
leptonic modes are allowed. If the splitting is smaller than the muon mass, only leptonic decays 
with electrons in the final state are allowed. We see, however, from table |IJ that this situation is 
not realized for our benchmark points. 



6.2. Charged tracks 

For very small mass splittings between the components of (f>, the (j)'^^ can travel over macroscopic 
distances before decaying. This is illustrated in fig. [6j which shows that mass splittings below 
Mn- are needed for the lifetime of (j)^ to become macroscopic. Note that the same relationship 
between AM and the lifetime cr applies also to the multiply charged components of (p. However, 
in our benchmark models, the smallest mass splitting and thus the largest cr is always the one 
corresponding to (j)^. 

If (p^ decays after travelling more than a few tens of centimeters, (p production can be potentially 
seen in searches for anomalous charged tracks in the inner detectors of ATLAS and CMS. Since 
in our benchmark scenarios, (p^ is part of all decay chains, such searches would be sensitive 
to the production of any charged component of (p. The cross section for this is very similar to the 
total 4> pair production cross section shown in fig. [sj Events with no charged (p (i.e. only (p^) are 
almost completely absent in our y = benchmark model, where cp^ does not couple to the Z, and 
they contribute only about 17% of the total cross section for the Y = 2 model. 

Searches for charged track signatures have been carried out by both ATLAS and CMS. We 
expect the best sensitivity to our benchmark models to come from future searches of the type 
presented by ATLAS in [83]. In this analysis, a high pT jet as well as more than 90 GeV of missing 
transverse energy are required in addition to the new charged particle. The latter is required to 
leave a signal only in the inner detector, making this search the most sensitive to particles with 
lifetimes on the order of several tens of centimeters. 
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Figure 6: Lifetime of the singly charged component of as a function of the mass sphtting AAf between (j) 
and The thick black curves show the physical lifetime, taking into account hadronic and leptonic decays. 
The colored curves show the inverse of the partial widths to pions (green), ei/g (red) and fiv^ (blue). The 
labels on the upper horizontal axis show the values of corresponding to the mass splittings indicated 
on the lower horizontal axis, neglecting electroweak corrections. 



The CMS search for long-lived charged particles [84J requires signals in the tracking detectors as 
well as the muon chambers, implying sensitivity only to particles with decay lengths of order 10 m. 
Similarly, the ATLAS search [85| requires signals in the inner detector and the electromagnetic 
calorimeter. Moreover, only particles with electric charges > 6e are constrained in this analysis. 
The searches from and [S^ are therefore not sensitive to our benchmark models or minor vari- 
ations therefore, except for an extremely fine-tuned corner of parameter space, where electroweak 
contributions to the mass splittings, eq. (jBJ), and those induced by nonzero X'^jj, eq. (5|, conspire to 
make one of the mass splittings extremely small. If we depart further from our bencnmark models, 
however, it is quite easy to obtain very long-lived charged particles. In particular, this is the case if 
the hypercharge is chosen such that the lightest component of (p is charged and decays only via 
higher-dimensional operators, for instance eq. ([s]). Then, its decay width is naturally very small. 
Note that in scenarios of this type, the long-lived charged particle should still decay on timescales 
<^ 1 minute to avoid perturbing big bang nucleosynthesis. 



6.3. Monophoton and monojet signatures 

Searches for a single jet or photon, accompanied by a significant amount of missing energy, 
have recently received a lot of attention because they are able to constrain the existence of new 
"invisible" particles in a relatively model-independent way |86l - ll03j . In the models discussed here, 
for instance, the components of the electroweak multiplet mediator 4> are very difficult to observe 
directly at the LHC, but their production is constrained by jet+MET and photon+MET searches. 

For the monojet signature, the relevant diagrams are pair productions of the multiplet together 
with a quark or gluon from initial state radiation. For our benchmark points the expected cross 
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Figure 7: The production cross section for a multiplet pair together with a monojet for our benchmark 
models. The cross sections (colored curves) are compared with a bound from the CMS monojet search [5S] 
(gray area), which requires the reconstructed MET to be above 350 GeV. Following [86 , we use a cut 
efficiency of 10% relative to a parton level Monte Carlo sample with parton level cut MET> 200 GeV. 
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Figure 8: The production cross section for a multiplet pair together with one photon for various multiplet 
quantum numbers. Based on the CMS 7 TeV monophoton search [ST], we require pt > 125 GeV and 
\r]\ < 1.5 for the photon, and we take the signal efficiency to be 30%. The shaded region is excluded by the 
CMS search. 
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Figure 9: Feynman diagram through which the new electroweak multiplet (j) contributes to Higgs annihilation 
into 77 and final states. 

section is still several orders of magnitude below the current LHC bound, as shown in fig. [7j 

The production cross section for a 4> pair together with a single photon is significantly enhanced 
compared to the monojet case because hard photons can be radiated not only from the initial 
state quarks, but also from the (j)'^^ in the final state, which couple strongly to photons. The 
monophoton cross sections for our benchmark points, as well as two additional cases with even 
larger SU{2) representations, is shown in fig. [sjand compared to the current limit from CMS |87j . 
We see that, in spite of the enhanced production cross section, the model is still not constrained 
by the current experimental searches. 

Finally, we have also considered the signature of (fxj)^ + h production through the operator 
^^(pH^H. The signature of this process — a Higgs boson plus a lot of missing energy — is identical 
to the one for associated Z + h production, in which the Z decays invisibly. For our = 5, 1^ = 
benchmark point, the cross section for (jxj)^ + h production at -y/s = 8 TeV is 0.3 fb, while for the 
= 5, = 2 benchmark point, it is 8 x 10~^ fb. Since the SM cross section for Z + hh production 
is ~ 400 fb, we do not expect to see any modification of the Higgs plus missing energy event rate 
in the foreseeable future. 



7. MODIFICATION OF HIGGS BOSON DECAYS 

In a model with multiple scalar fields, the presence of "Higgs portal" -type operators is quite 
natural; indeed, as these are dimension-four operators consistent with all of the other symmetries, 
they are difficult to forbid without ad-hoc assumptions. The presence of the operator (j)^(pH^H can 
significantly modify decays of the Higgs boson to gauge bosons, especially /i — )• 77 and h — )• Z7 
which arise in the SM only at loop level. 

We first consider the /i — )■ 77 decay width 



where v = 246 GeV is the vacuum expectation value of the SM Higgs and F"''^ a dimensionless 

777 
w ' 



amplitude. At one loop this amplitude receives SM contributions from the W boson loop, FT' 



and from fermion loops, Fj^ [66fl68[ riU4j . In our model there is an additional contribution, Fj'^ , 
from the scalar multiplet (p running in the loop (see fig. [9]), so that 

F^^ = F^^ + FJ' + . (27) 

Using analytic expressions from |67| [68] we have 

\F^\ = 1.25 X 10-=^ , = 2.75 x 10^^ , (28) 
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Figure 10: Ratio of partial decay width of the Higgs boson in our model to the standard model width for 
decay modes 77 (blue) and 7Z (red), as a function of the coupling A^//. The solid and dashed curves 
correspond to the stable and unstable benchmark points of table [T] respectively. Experimental results for 
best-fit signal strength in the channel h 77 (dashed horizontal line) are taken from |106l I107j with errors 
added in quadrature (gray band). 

for Mfj = 125 GeV, M\y = 80.4 GeV, and the top quark mass in the MS renormalization scheme 
Mf^^ = 160 GeV |78| . We neglect the contributions from fermions other than the top quark, which 
is sufficient to reproduce the full SM result for the partial width T{h — )• 77) |1U5| to within 2%. 
The new contribution Fj'^ is given by |66ft68j 

Z?77 



E ^ {^<^^ - [1 - Mm] , (29) 



where the sum runs over the components (j)s of the electroweak multiplet (p, and the loop function 
f{Ps) is given by eq. (|llj), and Ps = AM^^/M'^. For our benchmark models f3s > 1- Note that for 
y = 0, the term proportional to X'^^j vanishes. 

In the SM Fy^/ and Fj"' interfere destructively. Furthermore, for positive X(f,H and X'^jj- we have 

|^77|2 ^ (|^77| _ |^77| _ |^77|)2^ (3O) 

The ratio of the resulting partial width r(/i — t- 77) to the SM value as a function of X^pn is shown 



in fig. 10 We see that the decay /i — t- 77 can be substantially enhanced or suppressed, depending 
on the sign of the coupling constants X,f)H and X'^jj. Note that the other phenomenology discussed 
so far is to a large extent decoupled from the value of A^jy, so that 0(1) effects in /i — )• 77 are 
possible without affecting anything else. In particular, there is no clear prediction for the size 
of the deviation in /i — t- 77 based on observation of the Fermi gamma line, beyond the generic 
expectation that 0(1) deviation is expected for natural values of X^h- 

In a similar way, the related loop-induced decay h — )• 7Z is affected by the new multiplet (j)- 
The decay rate is given by 

r(/i^7Z) = ^^fl-^V|F^^|2, (31) 
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where the ampHtude F'^^ receives contributions from W loops, fermion loops, and loops, 

= + Ff + Ff . (32) 
Analytic expressions for these can be found, e.g., in refs. |67| I68j. In the SM we have 

= 1.41 X 10"^ (33) 
The new physics contribution is 



\F^\ = 2.63 X 10"^ 



V2^Qs((T3), - sl^Qs 



s 



Is 



+ 



2(/3, - 7,) 2(/3, - 7,) 



■[/(/3s)-/(7s)] + 



(/3s -7s)' 



[gm-g{ls)] 



(34) 



while (5s = 4M5 



I' 



Here, as in eq. (29), the sum runs over the components of 
and the loop functions f{(5s) and g(7s) have been defined in eq. (11) and eq. (12). Numerically, the 
new physics contribution is of the same order as the SM contribution for X^h ~ 0(1); see fig. 10 

Similar loop contributions correct the h — >• WW and h — )■ ZZ branching ratios. However, since 
these processes receive tree level SM contributions, the relative corrections from the new multiplet 
are small. They interfere destructively with the tree-level amplitude, leading to a slight reduction 
in the partial decay widths for h — )■ WW and h — )• ZZ. Numerical evaluation of the loop diagrams 
using FeynArts for our benchmark points yields corrections on the order of a few percent. 

Modifications of the other Higgs decay modes are negligible, since cj) has no direct coupling to 
any of the SM fermions. Invisible decays of the Higgs into </> or x are forbidden kinematically for 
the regions of parameter space we consider. 



8. CONCLUSION 

In conclusion, we have investigated the possible implications of gamma ray lines in astrophysical 
dark matter searches for the LHC. Motivated by the tenative hints for a line-like signal at ~ 130 GeV 
in Fermi-LAT data, we have focused on a class of secluded DM models, in which the DM particle x 
couples to Standard Model fermions and gauge bosons through loops of an intermediate particle (f). 
We have considered in detail the case where both x aiid 4' are vev-less scalars, but we expect our 
conclusions to be valid also in a more general context. We have moreover assumed that (j) belongs 
to a large representation (N > 3) of the weak SU (2) gauge group. 

Among the models proposed to explain the Fermi-LAT signal (if it stands up to further ex- 
perimental scrutiny), this scenario has several advantages: 1) For natural, untuned values of the 
coupling constants, it can explain relatively large {{av) ~ 10~^^-10~^^ cm^/sec) DM annihilation 
cross sections to 77 and/or 7Z final states. Both of these final states lead to monoenergetic features 
in the astrophysical gamma ray spectrum. The key is that some of the component fields of cj) carry 
several units of electric charge, which significantly enhances the DM coupling to photons. 2) All 
coupling constants are perturbative at experimentally relevant energy scales. 3) The model is well 
compatible with the observed DM relic abundance in the Universe. 4) DM annihilation to WW, 
ZZ, and fermion-antifermion final states is small. This is important because these final states are 
tightly constrained by searches for the broad excess they would induce in the Fermi-LAT data. 
5) DM-nucleon scattering cross sections are compatible with current constraints from direct DM 
searches, but are testable in future experiments. 

Turning to the LHC phenomenology, we found that the scenarios we consider have very char- 
acteristic signatures, but are still largely unconstrained at present. While the members of the 
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mediator multiplet cj) can be copiously produced due to their large electroweak couplings, they are 
difficult to observe because their cascade decays down to the lightest component field are typically 
very soft. The reason is that, unless there are large, isospin-dependent couplings to the Higgs, 
the mass splittings among them arise only from higher-order electroweak corrections. Thus, the 
energies of the decay products can easily be well below the trigger thresholds of ATLAS and CMS. 
The lightest component field, which we take to be the neutral one, (jP, in turn, is invisible due to 
its vanishing electric charge. In view of this, (j) production can contribute to final states with large 
missing transverse energy, for instance monophoton + MET. Here, the probability for radiating 
an extra hard photon in pp — )• (f)(j) is enhanced by the large electric charge of the component fields 
of <j). However, we have found that current monophoton + MET searches still fall two orders of 
magnitude short of constraining the most interesting regions of parameter space. 

A second place in which the multiplet can leave its footprint at the LHC is the Higgs sector. 
Higgs boson decays to 77 and 7Z receive extra contributions from diagrams involving (p loops, and 
these extra contributions can either enhance or suppress the corresponding Higgs branching ratios. 
The LHC data on /i — )• 77 thus already provides loose constraints on the (p-h couplings, as shown 
in fig. [10} 

Finally, if the mass splittings among the components of (j) are very small, some of them can be 
sufficiently long-lived to yield anomalous charged tracks that can be detected in future searches 
using the ATLAS and CMS inner detectors. Other signatures that might provide promising starting 
points for future work include photon -|- MET -|- X. With optimized cuts, these signatures could 
efficiently exploit the large (/>-7 couplings to improve the signal-to-background ratio. 
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Appendix A: SU(2) interactions 

In this appendix we give the normalization of the SU (2)^ generators used in the paper and then 
also write out explicitly the gauge interactions. For the generators of the representation N of 
the algebra sn(2), we fix the normalization by insisting that for any representation, the eigenvalues 
of (and thus the electric charge of the multiplet scalars) differ by integers. Thus, we have 
explicitly in the basis with diagonal 

iTff)mn = ^ ^ ^mm (^1) 

where m G [1,-^]- Since su{2) = so(3), the other generators of representation N can be obtained 
from the familiar angular-momentum ladder operators for spin j = (N — l)/2, 

(r+)mn = \/m{N - m)5m+l,n, {T^)mn = \/ {m - 1){N - m + (A2) 

and the relation = it iT^. It is easily verified that the three generators satisfy the 
defining relation 

[r^,r^] = ie'^'"^r^. (A3) 
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In the fundamental representation 2, these generators match on to the usual Pauli matrices, T2 
For the calculation of the T-parameter in section [5} we make use of the relations 



(^Af)L — (^Jv)ss'(^Af)s' 
1 



[s{N - s) + {s-l){N - s + l) 



(A4) 



and 



S + l)5s-l,s>) 



-^s{N - s)5,+iy + ^{s - 1){N - s + 



(A5) 



Using this normalization the Lagrangian eq. ^ is, in the gauge-field mass eigenstate basis, 



{W^{T+)ij + W;{T^)ij) + 5'^ -^Z^iT'^ - sin^ e^Q)ij + eA^Q 



cos ( 



+ 



4 



w^w^'^ {r+,r^}^. + 6'^ z^z^^L^^{T% - sin^ e^Qf + a^a^q\^ 



cos^ 



+2A^z^eg^(r3 - sin^^^Q)) + w^A^"^ ({t±,t3 },,^. + 2(r±),,y^) 

cos rj^ J y 2 

^W^Z^ J" ({T^, T^lii cos2 - 2{T^)i,Y^ sin2 6^) 



(A6) 



The A'^ji^ interaction term in eq. (jT| has an unusual form, and at first glance may not appear 
to be gauge invariant. We can demonstrate its invariance by performing an arbitrary SU(2) gauge 
transformation: 



b,d 



Expanding to first order in the parameter 9°", the bilinears transform as 



(A8) 



and similarly for H'^t°'H. The bilinear itself is not gauge invariant, but for the combination we 
find 



(A9) 



and the extra terms vanish by the total antisymmetry of the structure constants f°-^'^. 



Appendix B: Dark matter direct detection 



In this appendix we collect the analytical results for DM scattering on nuclei for our models 
where DM interacts with the visible sector through electroweak multiplets and the Higgs. The 
heavy fields, (p, W, Z and Higgs are integrated out and one matches onto the Effective Field 



22 



Theory (EFT) operator basis (13).^ These EFT operators induce spin independent DM-nucleon 



scattering, eq. (14). 



The only tree level contribution is due to a single Higgs exchange, giving 

Ml 

Cr = -Ax«4. (Bl) 

This contribution is absent in our benchmark points, where we set X^h = 0, but could in principle 
saturate the present DM-nucleon direct detection bounds if X^h ~ 3 • 10~^. For such small values 
of X^H the effects on annihilation cross section and early cosmology are atill very small, though. 
The Cp Wilson coefficients is first nonzero at 1 loop, where <j) fields run in the loop, giving 

with Qij) the charges of the (p field components, and the sum runs over all the components in the 
multiplet (for simplicity we have treated the masses of the (p components as degenerate) . For the 
5-plet benchmark model with = 0, we have thus Cp = 5aA^0/(127r) = 9 • 10~^, and for the 
5-plet with = 2 we find Cp = 5aX^^/{A7r) = 1 ■ W'^. 

The Cq Wilson coefficient also receives the 2-loop contributions with (p and W, Z, 7 in the 
loops. We calculate these contributions in the approximation where M,^ ^ ^w,z- In that case we 
can first integrate out the (p fields and match onto an EFT with the operator X^F^yF^^'^ from (13), 
as well as the operators 

^e// 3 %X'Z,.Z^^ + I^X^ Vi^^^ + ^x'W^yW-^'^ . (B3) 

' <t> 4> 

Here, we have defined Z^j/ = d^Zy — dyZ^ and W^^ = d^W^ — dyW^^. Note that electroweak 
symmetry is already broken in this EFT. The Wilson coefficients are 

247r [swcw) ~^ 



c^ = -^^TV(r+r-), (B5) 



Cz, = ^^Y.^Tl - sl,Q^)Q^ , (B6) 

with the sums again running over the components of (p. For the = 5, 1^ = benchmark model 
Cz = 3 X 10-^ Cvi/ = 8 X 10-3, = 2 X lO'^. For the N = f),Y^ = 2 benchmark model, the 
Wilson coefficients are = 2 x 10"^ Cvf = 4 x 10-^ Cz^ = 3 x 10"^. 

Integrating out W and Z, the one loop contributions with ZZ, WW and Z7 running in the 
loop give 

Cq = C^ + + Cp, (B7) 

with 

'^q = ^r^^ [ - '^Ql'w + '^QgSwTj^ + HTj^f + l2Qy^ log(M|//.2) [Tjj - Q^sl,]], (B8) 



Note that at this stage also the top quark should be integrated out. In what follows we will treat the top 
contributions to DM-nucleon scattering only very approximately. 
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-8^4;' ^ ^ 

^ « [ _ 5 + 6 log(M|/;x2)] [2Q,s'^ - T^] . (BIO) 

^ 47r swcw 

Appendix C: Calculation of the S- and T-parameters 

For the S'-parameter contribution from the scalar multiplct 0, wc need to compute the vacuum 
polarization amplitude n3y(0). At one loop, only a single type of Feynman diagram contributes, 
with an intermediate loop of scalar particles (if ^ had a vev, there would be additional contributions 
with internal gauge-boson propagators.) Labeling the external momentum as p and the loop 
momentum as k, the amplitude is given by the expression 

-T. j §f.w^iifT^,i-*s'^hA-2k-pr) ^ '_ ^ ,^ (-vy(2^+P)-). (CD 

Using the standard Feynman parameterization, we can shift the integration momentum and rewrite: 
TTM- . f\ f dH 4£/-r + (l-2x)Vp- 

The second term proportional to p^p*^ does not contribute to the transverse vacuum polarization 
(and thus to S), so we drop it. As for the first term, Lorentz invariance allows us to replace 
\i'^g^^ . Evaluating the momentum integral in dimensional regularization, we have 

nsY b') = E ^kss I' dxA,{l + E- log[A,//x2]), (C3) 

where = M| — x{l — x)p'^, E = 2/e — 7 + log(47r) — log(//^), and is the renormalization 
mass scale. It is clear at this point that if there is no mass splitting within the multiplet, then 
the amplitude is proportional to tr(T^) = 0, and there is no contribution to the S'-parameter. To 
convert to S, we need to take the "derivative" at p'^ = 0, i.e. 



S = -167rn'3y(0) = -167r lim 

p2->.0 



^{Usy{p')-nMo)) 



(C4) 



Making use of the identity 



log[{M^-xil-x)p^)/M^] _ -xjl-x) 



s 

we find that 

^ = E ^hs (2 + ^ + log(M,VM')) • (C6) 

s 

The first term vanishes due to the trace over T^, so the only contribution to 5" is due to the loga- 
rithm. Since the states of cf) come in pairs of equal but opposite Tff eigenvalues Is, the dependence 
on the scale /x cancels, and we are left with our final expression. 



sin ft... ma ft... ^-^ 



3 sin Qy, COS 



-fslog 



(C7) 
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The calculation of the T-parameter is somewhat more involved, since it depends on the cor- 
relation functions directly and not just their derivatives. This means that an additional diagram 
contributes to T, arising from the four-boson interaction (f)^ cpA^^a-^^'"" ■ Furthermore, the loop 
coming from three-boson vertices can now have two distinct species of in the loop, due to the 
off-diagonal structure of T^. We thus have 



(27r)4 A:2 _ + ie 



(C8) 
(C9) 



Once again carrying out the momentum integral in dimensional regularization, taking the transverse 
part, and evaluating at = 0, we find 



nL(o) 



'87r2 



/' dx[{xM^ + (1 - x)Mf,)(l + - log f^_Ml±0_^ 
Jo V M 



9 



Kaio) = E(^^)-^'[i + ^ - ^og{Ml/^?)]. 



(CIO) 
(Cll) 



Here {T^ is the square of the matrix element ss' of generator T^, not to be confused with the 
matrix element ss' of the squared generator {Tfj)'^. Making use of our explicit representation of 
the group generators, it can be verified that the 1/e divergence and scale dependence completely 
cancel in the difference nii(O) — 1133(0). The leftover contribution comes from nf]^(0), and is equal 
to 



nii(o) -n33(o) 



8^2 



4^ ' 2{M^-Ml) 



(C12) 



Making use of the identity eq. dASh, and the definition T = 47r/(s2m^)[nii(0) - 033(0)], we have 
finally 



-Q 



2sin4^^M^ 



E 



-^s{N - s)5s+i,s' + \{s - 1){N - s + 1)5,-1,,' 



log(M2/M2 



(C13) 



Appendix D: Tree level annihilation cross section 



If the coupling of DM to the Higgs is non-zero, it can annihilate at tree level annihilation to 
ZZ and fermion-antifermion final states which contributes to the astrophysical continuum 
photon flux. The cross sections for DM annihilation to gauge bosons are given by 
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where M^i/ [Mz) is the mass oiW (Z) boson. 

This gives {(TVrei)w+w- = ^ (2 • 10-^3 cm^/s), {(Jv^e\)zz = ^Ih ^ (1 • 10"^^ cm^/s) for 

our stable benchmark point {N = 5, = 0) and {o^v^■cl)^r+^,- = x (3 • 10~^^ cm^/s), 

{(Tv^c\) zz ~ '^H ^ (1 ■ 10^^^ cm^/s) for the unstable case {N = f),Y^ = 2). Using the continuum 
photon flux bound from the galactic center ~ 3 • 10~^^ cm^/s, this bounds X^^^h ^ 0.03. 
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